CONWAY'S NAPKIN PROBLEM 



ANDERS CLAESSON AND T. KYLE PETERSEN 

Abstract. The napkin problem was first posed by John H. Conway, and written up as 
a 'toughie' in "Mathematical Puzzles: A Connoisseur's Collection," by Peter Winkler. To 
paraphrase Winkler's book, there is a banquet dinner to be served at a mathematics con- 
ference. At a particular table, n men are to be seated around a circular table. There are n 
napkins, exactly one between each of the place settings. Being doubly cursed as both men 
and mathematicians, they are all assumed to be ignorant of table etiquette. The men come 
to sit at the table one at a time and in random order. When a guest sits down, he will prefer 
the left napkin with probability p and the right napkin with probability q = 1 — p. If there 
are napkins on both sides of the place setting, he will choose the napkin he prefers. If he 
finds only one napkin available, he will take that napkin (though it may not be the napkin 
he wants). The third possibility is that no napkin is available, and the unfortunate guest is 
faced with the prospect of going through dinner without any napkin! 

We think of the question of how many people don't get napkins as a statistic for signed 
permutations, where the permutation gives the order in which people sit and the sign tells 
us whether they initially reach left or right. We also keep track of the number of guests who 
get a napkin, but not the napkin they prefer. We find the generating function for the joint 
distribution of these statistics, and use it to answer questions like: What is the probability 
that every guest receives a napkin? How many guests do we expect to be without a napkin? 
How many guests are happy with the napkin they receive? 



1. Introduction 

The problem studied in this article first appeared in the book "Mathematical Puzzles: A 
Connoisseur's Collection," by Peter Winkler [Sj, and was inspired by a true story. Rather 
than recounting the problem and the story ourselves, we prefer to quote directly from "Math- 
ematical Puzzles": 

The Malicious Maitre D' 

At a mathematics conference banquet, 48 male mathematicians, none of them 
knowledgeable about table etiquette, find themselves assigned to a big circular table. 
On the table, between each pair of settings, is a coffee cup containing a cloth napkin. 
As each person is seated (by the maitre d'), he takes a napkin from his left or right; 
if both napkins are present, he chooses randomly (but the maitre d' doesn't get to 
see which one he chose). 

In what order should the seats be filled to maximize the expected number of 
mathematicians who don't get napkins? 

. . . This problem can be traced to a particular event. Princeton mathematician 
John H. Conway came to Bell Labs on March 30, 2001 to give a "General Research 
Colloquium." At lunchtime, [Winkler] found himself sitting between Conway and 
computer scientist Rob Pike (now of Google), and the napkins and coffee cups were 
as described in the puzzle. Conway asked how many diners would be without napkins 
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if they were seated in random order, and Pike said: "Here's an easier question — 
what's the worst order?" 

The problem of the Malicious Maitre D' is not horribly difficult; if you're having trouble 
finding a solution, you can see Winkler's book for a nice explanation. In this paper, it is 
Conway's problem that we focus on. Again, from [Sj: 

Napkins in a Random Setting 

Remember the conference banquet, where a bunch of mathematicians find them- 
selves assigned to a big circular table? Again, on the table, between each pair of 
settings, is a coffee cup containing a cloth napkin. As each person sits down, he takes 
a napkin from his left or right; if both napkins are present, he chooses randomly. 

This time there is no maitre d'; the seats are occupied in random order. If the 
number of mathematicians is large, what fraction of them (asymptotically) will end 
up without a napkin? 

Let p be the probability that a diner prefers the left napkin and q — 1 — p be the probability 
that a diner prefers the right napkin. For the case p — q — 1/2, Winkler's book gives two 
proofs of the answer: (2 — ^fe) 2 ~ .12339675. One is combinatorial, while the other, taken 
from a more general result due to Aidan Sudbury [3], is analytical. In fact Sudbury gives the 
expected proportion of diners without a napkin as (asymptotically) 

(l-pe^jl-qeP) 
pq 

(As an aside, we took an informal survey of 55 mathematicians and found about 69% would 
prefer the napkin on the left. According to Sudbury's result, we would thus expect about 
10.58 percent of the guests to get stuck without napkins.) 

In this paper, we use combinatorial methods to produce the generating function for the 
probability that at a table for n people, i of them have no napkin and j of them have a napkin, 
but not the napkin they prefer. This generating function allows for a thorough statistical 
analysis of the problem, including a new proof of 

^From our point of view, the number of people without a napkin is a statistic for signed 
permutations; just not one so well studied as inversions, descents, and such. We consider 
the order in which guests sit down at a place as a permutation of [n] = {1,2,..., n}, while 
their preference for the right or left napkin is given by a plus or minus sign. We label the 
places 1, 2, 3, . . . , n counter-clockwise, so that place % has place % — 1 on its left, place % + 1 on 
its right, and place n is to the left of place 1. With this convention, the signed permutation 
(2, —3, 4, —1) describes the following sequence of events at a table for four. The person sitting 
in place 4 sits down first and takes the napkin on his left. The person in place 1 sits next 
and takes the napkin on his right. The person at place 2 sits third and wants to take the 
napkin on his left, but since that napkin is already taken, he is forced to take the napkin on 
the right. Finally, poor person 3 sits last to find no napkin at his place. 

In Figure ^ we illustrate this sequence of events. There a "•" and a "A" represent a seat 
and a napkin, respectively. Moreover, a solid line symbolizes an arm taking a napkin and a 
broken line symbolizes an imaginary arm reaching for the preferred napkin. When looking 
at these diagrams, keep in mind that the table is circular: the guests at places 1 and 4 are 
neighbors. Each signed permutation thus corresponds to a particular set of diners who will 
have no napkin. In particular, for each permutation there is some number of diners without 
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Figure 1. A table for four. 

napkins, and we seek to determine what proportion of signed permutations leave a given 
number of diners napkinless. 

We remark that this problem makes for interesting mealtime conversation, and we have 
heard many suggestions for variations on the theme. What if instead of a math conference 
banquet with all men, it is a dinner party for couples, and the couples enter one at a time, lady 
sitting first? What if the dinner party is a mixer for singles, where now all the ladies enter 
first, sitting in alternating seats? 1 What if the guests don't mind looking farther afield than 
simply to their immediate left and right, say, reaching as many as two places over for a napkin? 
Exploring answers to these questions may lead to some interesting mathematics, but in this 
article we stay within the original question, where the guests are all male mathematicians, 
they enter in a totally random order, and they are too shy to reach beyond their immediate 
left or right to find a napkin. 

In section |2] we define the generating functions that we use throughout the paper. Section 
El explores the question of when everyone gets a napkin, and makes connections with com- 
binatorial objects called ordered bipartitions, due to Dominique Foata and Doron Zeilberger 
PP. This connection makes subsequent proofs much simpler, as in section HJ where we derive 
powerful identities involving our generating functions that ultimately lead to exact formulas. 
In section El we answer the original question of how many guests are expected to be without a 
napkin, as well as provide some other statistics of interest. Section |H1 gives a second derivation 
of the generating function for the expected number of guests without a napkin. 

2. Definitions 

Let p be the probability that a diner will reach to the left, and q = 1 — p be the probability 
of reaching right. We denote the set of signed permutations of [n] by C n . For any n £ C n , 
let \tt\_ and |7r| + be the number of negative and positive entries in tt, respectively. Also, let 
I 71 "! = I 71 "!- + = 7i be the length of tt. Let 0(77) be the number of guests without a napkin 
after every guest has been seated as described by n. Note that the number of guests without 

lr rhis variation is actually not very interesting mathematically, as it is not difficult to see the women will 
always get their choice of napkin, and that 1 /4 (or pq in general) of the men are expected to be without a 
napkin. 
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a napkin is equal to the number of napkins left on the table. Furthermore, let m(ir) be the 
number of people who get a napkin, but not their first choice. We say a guest is napkinless 
if he has no napkin, and a guest is frustrated if he gets a napkin, but not the napkin he 
originally wanted. 2 Otherwise, we say the guest is happy. Define the weight of n as 

tu(7r) :=p W -g W +x° W ?/ mW . 

For example, if tt = (2,-1,3,4), then \tt\- = 1, |7r|+ = 3, and \ir\ = 4. Also, the napkin 
between places 2 and 3 is unused (guest 4 is the unlucky one), and although person 1 gets a 
napkin, it was not the one he wanted. Thus, o(n) = 1, m(ir) = 1, and two of the guests are 
happy: w(n) = pq 3 xy. 

The generating function we are interested in is 

(2) C(p;x,y,z) := ^ S w ( n ">~ = Yl V T ( i >3,n)x t y 3 z n , 

n>l 7reC„ ' i,j>0 

n>l 

where pr(z, j, n) denotes the probability that at a table for n people, i of them are napkinless, 
and j of them are frustrated. 

The main result of this paper, Theorem|2l is an exact formula for C(p; x, y, z). Our approach 
is to first "straighten" the table. 

Suppose that instead of a circular table with n places and n napkins, we look at a straight 
table with n places and n + 1 napkins, so that each place has a napkin on both its left and 
right. If we know all that can possibly happen in this situation, then in order to determine 
the circular case, we just consider that the last person to enter the room sits "between" the 
first and last person on the linear table. Let us make this connection more precise. 

Let M n be the set of all signed permutations of [n] that result in taking neither napkins 
from the ends of the table. Similarly, C n (resp. lZ n ) denotes the set of all signed permutations 
of [n] that result in the left end napkin being taken but not the right (resp. right but not 
left), and B n denotes those signed permutations that result in both end napkins being taken. 
We note that C n = N n U C n U TZ n U B n . 

For the straight table we define the weight of tt as 

where o(ir) is still the number of people without a napkin at the (now linear) table, so 
that 0(71") + 1 is the number of napkins left on the table. Moreover, we define the following 



When discussing this problem, French mathematician Sylvie Corteel argued that if the mathematicians 
were French, they would never take the "incorrect" napkin. If the napkin they wanted was not there, they 
would simply cross their arms and refuse to eat. But of course, that is a different problem (how does the 
number of napkinless guests change as the proportion of French diners changes?). 
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generating functions for the linear table: 

n>0 7r£C„ 
n>0 ir£j\f n 

:= EE ; 

n>l 7re£ n 

:= EE ; 

n>l ir£lZ n 

:= EE ffi ^^ 

n>l 7rSB„ 

By construction, we have 

(3) 5(p; x, y, z) = N(p; x, y, z) + L(p; x, y, z) + R(p; x, y, z) + B(p; x, y, z). 

Further, we can observe that by symmetry 

C(p; x, y, z) = C(q; x, y, z), S(p; x, y, z) = S(q; x, y, z), 
N(p; x, y, z) = N(q; x, y, z), B(p; x, y, z) = B(q; x, y, z), 

and 

(4) R(p;x,y,z) = L(q;x,y,z). 

/From now on we will usually suppress the p and write C(x,y,z) instead of C(p; x,y, z), 
S(x, y, z) instead of S(p; x, y, z), etc. 

Now let us recast the generating function C(x,y,z) in terms of the generating functions 
for the linear table. Everything that has happened before the last person sits down at a table 
for n people can be considered the result of a signed permutation of [n — 1] playing out on 
a linear table. If the last person sits at place n, this is obvious, but in general the labeling 
of the seats is unimportant, i.e., we can always cyclically permute the labels on the seats 
without changing the weight of the permutation. 

In what follows, n = ar means that n is the concatenation of a with r. For instance, 
(3, -l)(-6)(2, -4, 5) = (3, -1, -6, 2, -4, 5). 

Let us consider a circular table for n guests and assume that 7r = o"(±n)r is a member of 
C n . Let 7r' = to. If the last person walks in and has both napkins available, then ir' must 
have resulted in leaving both end napkins on a linear table, that is, ir' G J\f n -i- Whether the 
last guest prefers the left napkin or the right napkin, he will get his choice. So, the weight of 
7r is: 

, . [ pwsin') if person n prefers left, 
wyK) — < 

\qws(n) if person n prefers right. 

If person n walks in to find only the left napkin available, n' G C n -±, and the last person will 
take that napkin regardless of preference, getting the one he wants with probability p, and 



S(p;x,y,z 
N(p;x,y,z 
L(p;x,y,z 
R(p;x,y,z 
B(p;x,y,z 
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Figure 2. C = z(N + (p + qy)L + (py + q)R + xB). 



getting frustrated with probability q: 

, . \pws(tc') if person n prefers left, 
| qyws{^') if person n prefers right. 

Similarly, if person n walks in to find only the right napkin available then 



w{ti) 



pyws(ir') if person n prefers left, 
qws(Tr') if person n prefers right. 



Finally, guest n can walk in to find both napkins already taken. In this case, we know that 
7r' G £> n _i and that guest n will be one of the napkinless guests: 



{pxws(^') if person n prefers left, 
qxws(ft') if person n prefers right. 



Figure El illustrates these, altogether 8, possibilities in the special case when the last person to 
sits at place n. All other cases can be reduced to this case by cyclically shifting the picture. 
In terms of generating functions we have showed that 

(5) C(x, y, z) = z(N(x, y, z) + (p + qy)L(x, y, z) + {py + q)R(x, y, z) + xB(x, y, z)). 
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3. A WARM-UP PROBLEM: WHEN DOES EVERYONE HAVE A NAPKIN? 

A natural question to ask (perhaps easier than the general question) is: what is the prob- 
ability, pr(0, l,n) = pr(0,n), that at a circular table for n people, every guest has a napkin? 
The generating function for this probability is C(0,y, z). By equation (JSJ), we have 

(7(0, y, z) = z(N(0, y, z) + (p + qy)L(0, y, z) + (py + q)R(0, y, z)) , 

but since (on a straight table with at least one person) there is always at least one person 
without a napkin if neither end napkin is taken, N(0,y,z) = 1, and the above equation 
reduces to 

(6) C(0, y,z) = z(l + (p+ qy)L(0, y, z) + (py + q)R(0, y, z)) . 

This makes intuitive sense because if everyone is to have a napkin, they all need to take the 
left napkin, or all need to take the right napkin. Therefore we turn our attention to L(0, y, z) 
(since R(0, y, z) follows by swapping p and q). 

For n a positive integer, let us consider a straight table for n guests and let tx G C n . Further, 
let 7i = <r(±n)r, where a and r are the signed permutations consisting of the letters to the 
left and to the right of n, respectively. 

If the last person to enter sits down at the rightmost seat (i.e., r is empty) then, since 
7r G C n , he necessarily prefers the left napkin. So, 

w s (n) = w s (cr)p. 

Assume that the last person to enter sits down with at least one guest to his right (i.e., r 
is not empty). Then he will be happy with probability p and frustrated with probability q: 

, , I wsMpwsir) if person n prefers left, 
\ws{<J)qyws{T) if person n prefers right. 

Thus the function L(0,y,z) satisfies the following differential equation: 

(7) ^L(0, y, z) = (L(0, y, z) + l) (p + (p + qy)L(0, y, z)). 

Se Figure El for an illustration. 

Let L n = [z n }L(l/2; 0, 1, z), where [z n ]F(z) denotes the coefficient of z n in F(z). Using flj|) 
we compute the L n for some small values: L\ = 1, L 2 = 3, L 3 = 13, L 4 = 75, L 5 = 541, 
and plug them into Sloane's Encyclopedia 0. Luckily, we get a hit with sequence 
A000670! These numbers happen to be fairly well known as the "ordered Bell numbers," or 
the number of ordered set partitions of [n]. We shall now give a bijection between the ordered 
set partitions of [n] and permutations for which everyone takes the left napkin. This bijection 
will also give the general solution to (J7|). 

3.1. Ordered set partitions. We now give a bijection between ordered set partitions of [n] 
and signed permutations that correspond to everyone at a linear table taking the napkin on 
the left. An ordered set partition, a, of [n] is a word 

a = B\B 2 ■ ■ ■ E>k 

where the "blocks" B\, B 2 , Bk are subsets of [n] such that {Bi, B 2 , . . . , -£>&} is a set 
partition of [n]. By convention we require that the elements of each block are written in 
decreasing order. 
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We describe the bijection with an example. Let 

« = {5,2}{6}{7,4,1}{3}. 

First, we give a minus sign to the least element of each block, then we remove the braces, to 
obtain 

n= (5,-2,-6,7,4,-1,-3), 

a permutation corresponding to a situation where everyone takes the napkin on the left. We 
can see that we indeed have produced a permutation in which everyone takes the napkin on 
the left by thinking of n as a diagram for the entry times and napkin preferences of everyone 
at the table. With 7r as above, we see that the person seated second from the right entered 
first and took the napkin on the left. The person seated third from the right wanted to take 
that napkin, but entered fourth, and so was forced to take the napkin on his left. Because 
we required that the elements of a block are in decreasing order, anytime someone has a 
preference for the napkin on the right, they find that it was taken by the person to their 
right. 

Clearly this process is reversible. Given a permutation where everyone takes the napkin 
on the left, it must have at least one minus sign. In particular, it must have a minus sign on 
1, since the first person to enter must take the napkin on his left. Anybody with a plus sign 
immediately to the left of a person with a minus sign must enter after that person. And if 
two or more people with plus signs are sitting (consecutively) to the left of a person with a 
minus sign, they must arrive in order of closeness to the minus sign; the closest to the minus 
sign first, followed by the second closest to the minus sign, and so on. This gives us the 
block structure of the partition on [n]. Reading tt from left to right, we separate the blocks 
by putting walls immediately to the right of any number appearing with a minus sign. For 
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example, if 

tt= (7, -5,-6,4,-1,3,-2), 
we convert this permutation into 

« = {7,5}{6}{4,1}{3,2}. 

Now we can build the generating function for L(0, y,z) by purely combinatorial means. 
Taking an approach from pQ (more on that paper in a bit), let 

pz pqyz i p(qy) n - 1 z n p(e qyz - 1) 

H(p; y, z):= — + — — + ■■■ + + 



1! 2! n\ qy 

Then H(y, z) := H(p; y, z) is the generating function for single blocks, {n, . . . , 2, 1}, since for 
every block, only the person corresponding to the least number gets the napkin he wants, 
leaving the other n — 1 people frustrated. For example with n = 5: 

/ X X X X X v 

w V V V \ 7 \ = p(w)* 

V qy qy p / 

Therefore 

(1 - H(y, z))- 1 = 1 + H(y, z) + (tf(y, z)f + ■■■ 
is the generating function for ordered sequences of blocks, and thus 

p(e qyz - 1) 



L(0,y,z) = (l-H(y,z))- L -1 



qy — p(e qyz — 1) 

One can easily check that this expression indeed satisfies the differential equation (JJJ). Due 
to equations (jlj) and (jUJ) we are now in a position to give the generating function for the 
probability that everyone gets a napkin on a circular table: 

C(0, 2/, z) = *(l + (p + gy)£(jo; 0, y, 2) + (py + g)L(g; 0, y, z)) 
(p + qy)p(e qyz - 1) (py + g)g(e py2 - 1) 



z 1 + 



92/ — p{e qyz — 1) PZ/ — q(e p y z — 1) 



3.2. Ordered bipartitions. Now we will generalize the correspondence just described. The 
paper of Foata and Zeilberger [T] introduces objects called ordered bipartitions, which are 
easiest to think of as ordered set partitions with some of the subsets underlined. A compatible 
bipartition is an ordered bipartition where all the underlined subsets are on the right, e.g., 

q = {5,2}{6} {!,4,7} {3}, 

where we adopt the convention that underlined subsets have their elements written in ascend- 
ing order. The bijection works as follows. For every non-underlined group in a, we perform 
the same operation as above, while for every underlined group we perform the "opposite" 
operation. Specifically, we put minus signs on all but the least element before removing the 
braces. This produces 

tt = (5, -2, -6, 1,-4, -7, 3), 

a permutation where everyone on a linear table receives a napkin. All we really need to 
observe is that the underlined groups correspond to the part of the table where people all 
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take napkins on their right, while the non-underlined groups all take napkins on the left. 
Thus, 

5(0, y, z) = (L(0, y, z) + 1) (R(0, y,z) + l) 

/ 8 n = vm?_ 

" (q(ePy z - 1) - py) (p(e<^ - 1) - qy) 

We can take this more general correspondence and see now that the permutations for 
which everyone takes the napkin on the left correspond to the ordered bipartitions with 
no underlined subsets, the ordered bipartitions with all subsets underlined correspond to 
permutations where everyone takes the right napkin, and the compatible bipartitions with at 
least one underlined subset and one non-underlined subset correspond to the permutations 
where everyone gets a napkin and both end napkins are taken. This last observation gives 

B(0,y,z) = L(0,y,z)R(0,y,z). 

4. Ordered bipartitions and generating functions 

Using ordered bipartitions (not simply the compatible ones), we can encode more than 
just those permutations where everyone gets a napkin. Using the algorithm below, we can 
encode the set of all signed permutations. Let ips be this map, where the subscript S reminds 
us that we are dealing with the straight table. The image of the injection ips will lead us 
to the main theorem of this section, which gives some wonderful relationships between the 
generating functions for the linear table. Given a signed permutation 7r of [n], we form its 
image fsi 71 ) as follows: 

(1) Find the least element, (ignoring signs), that is not already included in some 
subset. 

(2a) If ir(i) is positive, then underline the set including 7r(i), and set j — i + 1. 
While 1 7T (j? ) | > 1 7r(j — 1)| and ir(j) negative, 

add to the set containing and set j = j + 1. 
(2b) If ir(i) is negative, then set j — % — 1. 

While > + 1)| and 7r(j) positive, 

add ir(j) to the set containing ir(i), and set j = j — 1. 
(3) If every element is contained in a set, then delete all minus signs and quit. Else, go 
to (1). 

Clearly, no two permutations can be mapped to the same bipartition. We will demonstrate 
the injection with an example. Start with 

tt=(9,1,-3,2,5,6, -4,-7,8). 

The number 1 is the least element, and it is positive, so its set will be underlined. We start 
searching to the right of 1, looking for negative numbers with absolute value bigger than 1. 
We get 

9 {1,-3} 2 5 6 -4 -7 8 

Now 2 is the least element, and it is also positive, but there are no negative numbers imme- 
diately to its right, so it forms a singleton set, 

9 {1,-3} {2} 5 6-4-78 
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Now 4 is the least element not already in a set, and it is negative. So we start searching to 
the left of 4, looking for positive numbers with bigger absolute value. We get 

9 {1,-3} {2} 5 {6,-4} -7 8 

The next steps give, 

9 {l,-3}{2}{5} {6,-4} -7 8 

9 {l,-3} {2}{5}{6,-4}{-7} 8 

9 {l,-3} {2}{5}{6,-4}{-7}{8} 

{9} {l,-3} {2}{5}{6,-4}{-7}{8} 

Then we drop the minus signs to get, 

{9}{1,3}{2}{5}{6,4}{7}{8} 

so that we have an ordered bipartition of [n] where we write underlined sets in increasing 
order. 

We will now explain how the weight w(ti) of the signed permutation n can be read from 
the corresponding ordered bipartition a = ips(^)- To start with, the number of napkinless 
diners, o(n), is exactly the number of occurrences of an underlined set followed immediately 
by a non-underlined set. The frustrated diners are those who are not the least element in a 
block, less the people without any napkin: 

m(7r) = |a| — £(a) — o(tt), 

where |a| is the number of elements of the underlying set of a (same as \ir\ here), and £(a) 
is the number of blocks in a. Let a (resp. a) be the the bipartition formed from the blocks 
of a that are underlined (resp. non-underlined). For underlined blocks of a we have that 
the least element is positive in tt and the other elements are negative in it. Similarly, for 
non-underlined blocks of a, the least element is negative in 7r and the other elements are 
positive in it. Thus, 

I 71- !- = M - ^(«) + ^(«); 

|7r|+ = |a| — £(a) + £(a). 



Let (ps{C n ) be the image set of bipartitions corresponding to all signed permutations of [n]. 
By examining the algorithm describing ips, we see that the set (ps(C n ) consists of all ordered 
bipartitions that never contain the patterns 

■ ■ ■ a] {b} ■ ■ ■ or • • • {b} {a ■ ■ ■ , where a < b. 

To simplify notation in what follows, we write C for C(x, y,z), S for S(x, y, z), etc. Further, 
we write H for H(p; y, z) and E_ for H(q; y, z). The following theorem takes advantage of our 
combinatorial model. Its power is reflected in the subsequent results it implies. 

Theorem 1. We have the following formulas: 

(9) L + B = HS 

(10) B — HSH 
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Proof. If we add a non-underlined block of size r to the left of any bipartition in (ps{C n ), 
then clearly we get a bipartition in (fs(C n+r ) that corresponds to a permutation in C n+r U 
B n+r . Furthermore, this new block will not change the number of people without a napkin 
(occurrences of underlined blocks immediately to the left of non-underlined blocks), and the 
number of new people who get a napkin they don't want is exactly r — 1. Therefore, 

L + B = HS. 

Similarly, if we add a non-underlined block of size r to the left of a a bipartition in <ps{C n ), 
and an underlined block of size s to the right, then we get a bipartition in ips{C n+r+s ) that 
corresponds to a permutation in B n+r+s . Thus, 

B = HSH, 

which concludes the proof. □ 
Corollary 1. We have the following formulas: 
N = S(l-H)(l-H) 
L = SH(l — H) 

R = S(l-H)H 
B = SHH 

C = zS{\ + q(y - 1)H + p(y - l)H + (x - y)HH) 

Proof. From Theorem^ we immediately get the formulas for B and L. Plug these formulas 
into (J3J) and Q and the formulas for R and N follows. Finally, (0) yields the formula for 
C. □ 

What Corollary n tells us is that if we can find an explicit formula for S, then we will have 
explicit formulas for all the other generating functions. We will derive such an explicit formula 
shortly. First we need to introduce the notion of a cyclic bipartition. A cyclic bipartition is a 
bipartition for which one element is distinguished, and only the cyclic ordering of the blocks 
matters. As a convention, we put the block containing the distinguished element at the far 
right if that block is not underlined, or at the far left if it is underlined. In our notation we 
will enclose the distinguished element in parentheses. For example, 

{1,6} {8,2} {3,5} {9,(7),4} 

is a cyclic bipartition, which we also write 

7,4} {1,6} {8,2} {3,5} {9, 

so that the distinguished element is equivalently the first element in a bipartition where a 
block can "wrap around." Similarly to how we encoded signed permutations playing out on 
a straight table with ordered bipartitions, we can encode the case of the circular table with 
cyclic bipartitions. Here the distinguished element will correspond to a "distinguished guest," 
who sits in place number 1 at the circular banquet table. 

Let ip be the map encoding the circular table case. For any signed permutation 7r of [n], 
we form its image (p(ir) as follows (the only difference in our algorithm is that the searches 
are cyclic): 
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(1) Find the least element, (ignoring signs), that is not already included in some 
subset. 

(2a) If 7r(i) is positive, then underline the set including ir(i), and set j — i + 1 mod n. 
While > — 1)| and ir(j) negative, 

add 7r(j) to the set containing 7r(z), and set j = j + 1 mod n. 
(2b) If 7r(i) is negative, then set j = i — 1 mod n. 
While |vr(j)| > |7r(j + 1)| and ir(j) positive, 

add ir(j) to the set containing 7r(z), and set j = j — 1 mod n. 
(3) If every element is contained in a set, then delete all minus signs and quit. Else, go 
to(l). 

As an example, start with the permutation 

vr= (-7,1,-3,4,-2,5,-6). 

The steps of the encoding are 

-7 {1,-3} 4 -25-6 

-7 {l,-3} {4,-2}5 -6 
and this last step, which requires us to perform a cyclic search, 

-7}{l,-3} {4,-2} {5,-6, 
Now we can drop the signs to get 

7} {1,3} {4,2}{5,6, 

or 

{5,6,(7)}{1,3} {4,2}. 

Let (p(C n ) be the set of all cyclic bipartitions corresponding to signed permutations on a 
circular table. We will use this set of cyclic bipartitions along with ipsiCn) to obtain the 
following theorem regarding the derivative of S(x,y,z). 

Theorem 2. We have the following formula relating the generating function for the straight 
table and the generating function for the circular table: 

(11) z^-S{x, y, z) = C(x, y, z)S(x, y, z). 

dz 

Proof. To prove equation (fTTj) . it will suffice to equate coefficients and prove 

(12) nS n {x, y) = ^2\ n .) C k {x, y)S n „i{x, y), 

i=i \ l ' 



in which the polynomials 



S n (x,y) := ^ w s{n) and C n (x,y) := ^ 



W(TT J 



are the coefficients of z n /n\ in S(x,y,z) and C(x,y,z), respectively. We will prove (fT^j) 
bijectively. 

The left-hand side of equation ()12|) can be thought of as counting the weights of permu- 
tations corresponding to bipartitions in (p s (C n ) with a distinguished element, or a straight 
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table with a distinguished guest. We simply put parentheses around one of the n elements of 
the bipartition, as in 

{2,8} {1}{7,(4),3}{9} {5,6} . 

Given any bipartition with a distinguished element, we can form a pair (c, s), where c 
is a cyclic bipartition of a subset A C [n] (corresponding to one in ip(C n )), and s is an 
ordered bipartition of [n] \ A (corresponding to one in v?s (£«))• We simply split the table 
into two pieces just before or after the block containing the distinguished element. If the 
block containing the distinguished element is not underlined, we make the split just after 
that block. The above example yields the pair 

({2,8}{1}{7,(4),3},{9}{5,6}). 

If instead the block with the distinguished element is underlined, as in 

{2,8}{1} {3,(4),7} {9} {5,6} , 

we make the split before the underlined block. Now the right half becomes the circular table, 
and the left half is the straight table: 

( {3,(4),7} {9}{5,6},{2,8}{1}). 

By splitting the bipartition as we do, all of the guests retain their status as happy, frustrated, 
or napkinless. In other words, the product of the weights of the pair of tables equals the 
weight of original table. 

Now, given any pair (c, s), where c G f(Ci) and s G (ps(C n _i), we can form a bipartition 
in (ps(C n ) with a distinguished element. First, in any of (") ways, we choose a subset A = 
{ai < a 2 < • • • < a^} C [n], and replace k with a k in c. For s, we replace k with b k , where 
[n] \ A = {bi < b 2 < ■ ■ ■ < b n -i}. Now we concatenate the bipartitions. If the distinguished 
element of c is in an underlined block we put the straight table on the left: sc. If the 
distinguished element is not in an underlined block, we put the straight table on the right: 
cs. □ 

Now, thanks to Corollary ^ and Theorem |21 we have: 

3'=({x- y)HH + q(y - 1)H + p(y - l)H + l) S 2 . 

We can solve this differential equation to get the following exact formula: 

(13) S(x,y,z) = —p-, 

where the denominator D is: 

pq(y — x)e yz + (qx — pqy — q 2 y 2 )e pyz + (px — pqy — p 2 y 2 )e qyz 

(14) 

+ pq(y{y-l) +x(l-yz))+y - x. 

We can now obtain exact formulas for any of the other generating functions discussed here 
by plugging equation (fT3~j) into the formulas of Corollary [TJ In particular, we have our main 
result. 

Theorem 3. At a table for n people, the probability that i people are napkinless and j people 
are frustrated is given by the coefficient of x l y^z k in the following function: 

pqyzUx — y)e yz + (qy 2 + py — x)e pyz + (py 2 + qy — x)e qyz + x) 
(15) C{x, y, z) = , 
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where D is given in ()14j) . 

5. The expected number of napkinless guests and other statistics 

With the generating function from (fL3j) . we can in principle extract any sort of statistics 
related to the number of napkinless guests or frustrated guests. We will highlight a few 
statistical results that are of interest to us. In particular, we find the expected number 
of napkinless and frustrated guests, the variance for each of these distributions, and the 
covariance for their joint distribution. 3 

5.1. The expected number of napkinless guests. Recall the definition of the polynomial 
C n (x, y) = J2-K€C n M 77 ) from the proof of Theorem[21 Suppose we know exactly what C n (x, y) 
is for some n. If we want to obtain the expected number of people without a napkin, we want 
to compute the weighted average 

E n (o(7c)) := ^pW- g H+o(7r). 

In terms of the polynomial C n (x, 1) (we set y — 1 since we're not interested in frustrated 
diners at the moment), this just means that we differentiate with respect to x and set x — 1, 
since 

C' n (x,l)= P |7r| -g |7r|+ o(7r)x°W- 1 

ci(i,i)=Ef H i H+ °W- 

Therefore, we want to find the generating function for the numbers E n (o(Tr)) = C' n (l, 1), 



or 



E(z) :=^E n (o(7r))z" 

n>0 

d 



z(pq{2 — z)e z + (p 2 + pqz — l)e pz + (q 2 + pqz — l)e qz + l) 



x=l 



(16) ... ., 

pq(l — zy 

Theorem 4. The expected number of napkinless guests on a circular table for n people is 

try 

(17) E n (o(n)) = — (l -pexp n (q) - qexp n (p) +pqexp n (l)), 
where exp n (x) = X]fc=o x V^! ^ s ^ e truncated exponential function. 

Proof Let f(n) = nil — pexp n (q) — qexp n (p) + pq exp n (l)) / (pq) be the right hand side of 
equation (|T7jl . and let F(z) be the ordinary generating function for the numbers fin). Note 
that 

nf(n + l) = {n + l)f{n)+ 1 ~ ''" 



n - l)\ 



3 This task is daunting by hand, but luckily we have technology to help us. Computer software such 
as Maple, for example, is very helpful, both with solving the differential equation leading to l|13fl . and in 
obtaining residues. 
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which implies 



d_ 

dz 



or equivalently, 



[zF(z)] + z{e z - pe pz - qe qz ) 



I - : )]■'( : ) f I f - ) F(z) + z(e z - pe pz - qe qz 



z 



It is easy to check that E(z), given by formula ()16j) . satisfies this differential equation. □ 

Formula (|T7|) implies Sudbury's result (equation (JTJ)). 
Corollary 2. The expected value 0/0(71") satisfies 

n , , ,v n(l -pe 9 )(l - ge p ) / 1 

\n! 

In particular, when p = q = 1/2 we have 

E n (o(n)) =n(2-v^) 2 + oQ). 

So, the answer to Winkler's problem of napkins in a random setting is (2 — \fe) 2 ~ 
0.12339675. It was quite a bit of work for this answer, but of course our work pays off 
in being able to find the following statistics as well. 

5.2. Other statistics. In the rest of this section we present asymptotic estimates for further 
statistics regarding the napkin problem. Since the formulas are less messy, we will primarily 
restrict our attention to the p = q = 1/2 case, but our approach is general. 

It is straightforward to obtain asymptotic estimates for functions with a finite number of 
poles such as E(z). See, for example, chapter 5 of Herbert Wilf's book We use the same 
technique to obtain all the estimates given here. We briefly outline the approach. 

Suppose we have a power series f(z) = Y^ o, n z n , with a singularity at z = 1 of multiplicity 
m (and no other singularities). Then the Laurent expansion of / around 1 is: 

f(z) = J^ + ... + ^ + b + bl (l-z) + b 2 (l-zf + .... 

If we let 

/ \ \^ n b- m 6_i 

q[Z) = > C n Z = -, r h • V 7 r, 

(called the principal part of /), then the function h(z) = f(z) — g(z) is entire, and its coeffi- 
cients vanish very quickly. Thus for large n, the coefficient of z n in g(z) closely approximates 
the coefficient of z n in f(z), i.e., 

fm + n — 1\, /n + l\. , 

a n ~ c„ = I I o_ m H hi I o_ 2 + o_i. 

As an illustration, we apply this technique to E(z) to, again, derive Sudbury's result. 
We first expand E(z) as a series in it = 1 — z to get 



E(u) = ^ b n u n 

n>-2 

(1 - u) (pq(l + n)e 1 " u + (p 2 + pq(l - u) - l)^ 1- ^ + (q 2 + pq(l - u) - l)e 9(1 " u) + l) 



pqu 2 
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We let E{u) := u 2 E(u), then set u = to obtain 

-r/^ , PQ e + (P 2 +pq - l)e p + (q 2 +pq - l)e g + 1 
Liu) = 0-2 = • 

pq 

Next, we differentiate E before setting u — 0. This gives 

. pg(e + e p + e 9 ) + (l+p)(p 2 +pg-l)e p + (l + g)(g 2 +pg-l)e 9 -l 
-C/ (0) = 0_i = . 

pq 

Now we can obtain our estimate: 

a n ~ (n + 1)6_2 + 6_i = — — — — . 

pq 

Now, if we want to get the variance of the distribution of napkins on the table, we need to 
compute the sums of the squares of the number of napkins left on the table, since 

Var n (o(7r)) = £„(o(tt) 2 ) - E n {p{ir))\ 

and we already have E n (o(%)). We get 

E n (o(7i) 2 ) = p w -? w+ o(7r) a 

= c:(i,i)+^(i,i). 

Therefore the generating function for the second moment is 

E n {o{K) 2 )z n = ^ [C(x, 1, z)] ^ + E(z), 



n>0 



2 



d 

and all we need to find is -—r 

dx z 



C(x,l,z) 



x=l 



Computing, we find: 



2z(2 - e z ' 2 ) (e z ' 2 (l - zf + (2 - e z ' 2 ) {(e z ' 2 - 1) 2 (3 - z) - e z + 2z)) 
(7^(1/2; 1, 1, z) = ^ ^— L 

Using the same method, we find that the variance is asymptotically 

n(l — pe q )(l — qe p ) (l — {p 2 — pq)e q — (q 2 — pq)e p — pq(e + 1)) 



Var n (o(7r)) 



p 2 q 2 



or n(3 - e)(2 - y/e) 2 w n(.0347631) for the p = g = 1/2 case. 

To get the expectation and variance for the number of frustrated guests, we follow the same 
procedure, except now we differentiate C(l,y, z) with respect to y. The covariance is 

m(7r)), 

so we differentiate with respect to x, then with respect to y to get the generating function 
for £'(o(7r)m(7r)), the only piece we don't know. The statistics are summarized in Tabled 
Notice in particular that if the expected number of napkinless guests is n(2 — y/e) 2 , and the 
expected number of frustrated guests is n(6y/e — e — 7), then the expected number of happy 
guests is n(4 — 2y/e) ~ n(. 702557). Seventy percent of the guests are happy! 
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X: 


o(ir) (Napkinless) 


i7i(tt) (Frustrated) 


E(X) 


n(2 - v^) 2 w n(. 12339675) 


n{6^-e-7) « n(. 174046) 


Var(X) 


n(3 - e)(2 - v^) 2 « n(. 0347631) 


^(6^ - e 2 - e - 38^ + 46) 
« n(. 13138819) 


Cov(X, Y) 


n ( - (2 - y/e){V<? - 3e - 5^ + 12)) « n(-. 029239461) 



Table 1. Statistics for napkinless and frustrated guests with p = q = 1/2. 



6. Another proof for the expected number of napkinless guests 

Upon reviewing a draft of this paper, Ira Gessel pointed out that the generating function 
for the expected number of napkinless guests satisfies the following differential equation (with 
P = q = l/2): 



where N(l, 1, z) is the generating function for the proportion of signed permutations for which 
neither napkin on a straight table is taken. We will think of such permutations in terms of 
their image under (p$'- specifically, those ordered bipartitions in for which the leftmost block 
is underlined and the rightmost block is not underlined. Gessel suggested that there may be 
a simple combinatorial explanation for (|18|). and there is. 

First, if all we want is the expected number of napkinless guests on the circular table, then 
because of the symmetry of the table, we have E n (o(ir)) = np 1 , where p' is the probability 
that any particular guest (say guest 1) has no napkin. Upon equating coefficients, equation 
(JTHJ) claims that \M n \ = \C*\, where C\ is the set of all signed permutations for which guest 1 
gets no napkin on a circular table. 

We can give a bijection between the ordered bipartitions in tps{N n ) and the cyclic bi- 
partitions in <p(Cn) as follows. Given an ordered bipartition for which the leftmost block is 
underlined and the rightmost block is not underlined, there must be a guest in the middle 
of the table who is napkinless. Make the leftmost such guest distinguished, and cyclically 
permute the blocks of the bipartition until the block with the distinguished guest is first or 
last, depending on whether that block is underlined. For example, 

{8}{(9),6,1}{2} {3,5} {7,4} 

is an ordered bipartition in ips(J\f n ) with the leftmost napkinless guest highlighted (guest 2 
in this case). We cyclically permute the blocks to get 



E(z) 



a 

dz 



1 - z 



a 

dz 



N{l,l,z) 



{2}{3,5}{7,4}{8}{(9),6,1}, 
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which is a cyclic bipartition where guest 1 gets no napkin on a circular table. The inverse 
of this bijection is given by simply cyclically permuting the blocks until we have an ordered 
bipartition for which the distinguished guest is the leftmost guest without a napkin, the 
leftmost block is underlined, and the rightmost block is not. It is straightforward to check 
that we can always achieve this state. 
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